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The formation of bound states of holes in an antiferromagnetic spin-1/2 back- 
ground is studied using numerical techniques applied to the t — J Hamiltonian on 
clusters with up to 26 sites. An analysis of the binding energy as a function of clus- 
ter size suggests that a two hole bound state is formed for couplings larger than a 
"critical" value J/t] c . The symmetry of the bound state is d x 2_ y 2. We also observed 
that its "quasiparticle" weight Z2h (defined in the text), is finite for all values of the 
coupling J/t. Thus, in the region J/t > J/t] c the bound state of two holes behaves 
like a quasiparticle with charge Q = 2e, spin 5 = 0, and d x 2_ y 2 internal symmetry. 
The relation with recent ideas that have suggested the possibility of d-wave pairing 
in the high temperature cuprate superconductors is briefly discussed. 
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The study of high temperature superconductors continues attracting considerable at- 
tention. Recently, a novel theoretical scenario has been proposed (and supported by self- 
consistent calculations, where the symmetry of the superconducting condensate is d x 2_ y 2, 
instead of the more standard s-wave of the BCS theory. On the experimental side, the Lon- 
don penetration depth, A(T), has been measured at low temperatures in YBa 2 Cu30 6 .9 5 . 
A linear temperature dependence compatible with a d-wave superconducting state was ob- 
served. Angular resolved photoemission experiments in Bi2Sr2CaCu208+<5 have found an 
anisotropic superconducting gap also compatible with d x 2_ y 2 superconductivity. [[| Before 
these recent developments, the presence of an attractive interaction in the d x 2_ y 2 channel 
appeared frequently in the theoretical analysis of holes in antiferromagnetic backgrounds. 
It has been shown that the interchange of magnons between carriers naturally leads to an 
attraction that is dominated by the d-wave channel. |3J Numerical studies have also consis- 
tently suggested that two holes on finite clusters form a bound state with d x 2_ y 2 symmetry 
in a staggered spin background. However, it is not clear if these results are an artifact 
of the approximations made in their calculation. Since in the cuprates superconductivity 
appears in the presence of antiferromagnetic correlations, a connection between these results 
and real materials may exist, and thus it is important to understand if the above mentioned 
d-wave attraction indeed exists in realistic models of correlated electrons. 

The purpose of this paper is to report on a detailed numerical study of the formation of 
hole bound states in an antiferromagnetic background represented by the t — J model. A 
study of their properties is presented for several cluster sizes in the subspace of zero and two 
holes, which allow us to analyse finite size effects. Dynamical properties of the two holes 
bound state are discussed, including the "quasiparticle" weight, Z 2 h, obtained by creating a 
pair of holes with the appropriate rotational and translational symmetry over a fluctuating 
spin-1/2 antiferromagnetic background. Our results for different cluster sizes confirm that 
two hole carriers on such a spin background tend to form a bound state in the d-wave spin- 
singlet channel for couplings J/t larger than a critical (finite) value J/t] c . In this regime, 
individual holes are unstable towards pair formation, and the bound state behaves like a 
dressed quasiparticle with charge Q=2e, spin and internal d-wave symmetry. 



In addition, the present paper also addresses an important issue that has been highly 
controversial in the context of theories of strongly correlated electrons. Does the quasiparticle 
weight of holes injected in an antiferromagnetic background vanish? In other words, do 
holes behave like quasiparticles? While these questions can be experimentally settled by 
photoemission techniques in single crystals of cuprate superconductors, the analysis of the 
results is difficult due to the presence of a considerable background in the signal. Thus 
a definite answer has not been experimentally given. || In some theoretical approaches, 
like the RVB scenario, the dressed holes ("holons") carry charge e but zero spin, and the 
missing spin 1/2 is carried by a charge neutral excitation called spinon. [jllj These ideas 
are supported by calculations in one dimensional models where the separation of spin and 
charge indeed takes place. In this toy model, the wave-function renormalization of one hole, 
Zih, vanishes in the bulk limit at the Fermi surface, and thus the Fermi liquid fixed point 
does not exist in one dimension (ID). Jl0| Although, there is no reason to assume that one 
and two dimensions are qualitatively similar, nevertheless it is important to explore this 
exotic scenario. The marginal Fermi liquid theory also interprets the photoemission results 
in terms of a vanishing Z^. |TT| On the other hand, novel but more conservative ideas, 
like Schrieffer's spin-bag approach, [TI| predicts the reduction of Z lh due to the presence 



of strong correlations to values considerable smaller than those for a weakly interacting 



system, but remaining finite in the region of physical interest. [13| Such a reduction in 
affects the interpretation of the experimental data, but does not change the basic idea of the 
pairing theory where dressed quasiparticles form bound states at low temperatures due to the 
interchange of some suitable excitation. Thus far, numerical studies of Zxh have confirmed 
these ideas, namely that Zih is small but finite in the interesting region of parameter space. 



[p~3|j In particular, a recent finite size study on clusters of up to 26 sites has provided strong 



evidence in favor of this result. ]14]] However, due to the availability of only a discrete set of 
momenta, the currently available lattices larger than a 4 x 4 cluster do not allow the study 
of holes right at the Fermi surface |IJ|. Here, we address this problem by a study of the 



wave-function renormalization in the subspace of two holes, since in this case the ground 
state of two holes always belongs to the zero momentum subspace which is contained in 



all clusters we have analyzed in the present paper. In agreement with the previous results 
obtained for one hole, in the present study we found that Z 2 h is finite for all values of the 
ratio J/t different from zero. [|1^ 

The t — J model is defined by the Hamiltonian, 

H = J £(S,.Sj - I ninj ) - t £ (cl iS c jiS + h.c.), (0.1) 

(y> <ij),s 

where cf s denote hole operators; nj = + n^; and clusters of N sites with periodic bound- 
ary conditions are considered. The rest of the notation is standard. The calculations have 
been carried out on square clusters satisfying N = n 2 + m 2 (where n, m are integers, and N 
is the number of sites) with N = 16, 18, 20 and 26 sites. These clusters are commonly studied 



to search for a smooth extrapolation of the results to the bulk limit. [17] The algorithm used 
was the standard Lanczos method, using translational symmetry to reduce the size of the 
(sparse) Hamiltonian matrix, as well as rotations in 7r/2 and spin inversion. For the N = 26 
cluster the Hamiltonian block corresponding to the ground state (GS) symmetry sector (Bi) 
with the largest size has 4,229,236 states, which can be handled only by supercomputers like 
the Cray-2. |Jj 

Intuitively, it is clear that a bound state of two holes will be formed at large values of 
J/t. The reason is that each hole individually "breaks" four antiferromagnetic (AF) links, 
which costs an energy of the order of the exchange coupling. Two holes minimize the lost 
energy by sharing a common link, and thus reducing the number of broken AF links from 
eight to seven. When the coupling J/t is reduced to more realistic values, this attraction 
may survive till some critical coupling is reached. A smooth behavior of the ground state 
energy of two holes (measured with respect to the zero hole energy) is observed both as a 
function of the coupling and of the cluster size, |RJ suggesting that our results for the energy 
are close to the bulk limit. Having calculated the exact ground state by Lanczos methods, 
the symmetry of the two holes ground state can be easily analyzed. For all clusters studied 
here it has been found that the ground state belongs to the Bi irreducible representation of 
the C 4 point group of the square lattice, which is equivalent to the d-wave symmetry. fi~9 | 



Note that, since the symmetry group of the clusters does not, in general, include reflexions 



we cannot exclude a partial mixing between d x 2_ y 2 and d xy symmetries (except for 16 and 
18 sites). 

In Fig. la, the average distance between the two holes, obtained from a study of hole- hole 
correlations in the exact ground state wave functions, is plotted as a function of J/t for 
different cluster sizes. It is observed that for coupling J/t = 0.5 or larger, the results seem 
to have converged to a finite number, indicating the presence of hole binding. On the other 
hand, the results at J/t = 0.2 show that the average hole distance grows appreciably as the 
lattice size is increased, and thus binding may not occur in this regime. We believe that in 
the intermediate region a critical coupling exists where binding of holes starts. 

In Fig. lb the binding energy of two holes is presented. This quantity is defined as 
A# = &2 — 2ei, where e n = E n — E , and E n is the ground state energy of the t — J model 
in the subspace of n holes. j2(J The single GS energies e\ have been calculated elsewhere 



on clusters up to 26 sites. [[HJ If < in the bulk limit, a bound state of two holes 
is formed. is an intensive quantity and thus it is more severely affected by finite size 
effects than the (extensive) ground state energy. In addition, the energy of one hole enters in 
the definition of A B and, as discussed before, this quantity carries an additional systematic 
error due to the absence of momentum k = (tt/2,it/2) in the discrete set of momenta of the 
clusters with N = 18, 20 and 26 sites, used in the present study. In spite of this problem, 
qualitative information can be obtained from Fig. lb with some confidence. The "critical" 
coupling, J/t] c where two holes reduce their energy by forming a bound state, slowly grows 
with increasing lattice size suggesting that it may converge to a finite value in the bulk 
limit. In Fig. lb, recent Green Function Monte Carlo results on 8 x 8 clusters are also shown. 



21|] With this approach, supplemented by the use of appropriate variational guiding states 
obtained from the analysis of smaller clusters, it has been possible to study couplings J/t 
as small as 0.4. The dashed line in the figure shows an educated extrapolation suggesting 
that binding starts at J/t] c ~ 0.3, in qualitative agreement with our exact results on smaller 
lattices. |22| Note also that calculations using a recently developed "truncation" method have 
provided evidence that the t — J z model, i.e. a model where transverse spin fluctuations are 
switched off, has also a finite critical coupling beyond which two holes form a bound state. 



In this model, J z /t] c ~ 0.18, which is in qualitative agreement with our results for the t — J 
model which suggest a larger critical coupling, since in the absence of spin fluctuations a 
stronger tendency to pairing would be expected. |f23[ 

In Fig. 2a, the wave-function renormalization Z 2 h is shown as a function of 1/N for several 
coupling constants. If two holes form a bound state then the analysis of Zih becomes irrel- 
evant, since isolated holes will become unstable towards pair formation. In other words, in 
simulations carried out in the grand canonical ensemble, where a chemical potential \x selects 
the fermionic density, the state of one hole is never stable in the region of pair formation. 
In such a regime the elementary charge carrier will be the two holes bound state, and thus 
it is necessary to analyze Z 2h , defined as, 

Z 2h = Wgsl aWgsn , (0.2) 

^|a«a£|v#> 

where \?pgs) is the ground state in the subspace of n-holes which can be obtained using 
exact diagonalization methods. The operator that destroys a pair of holes is defined as 
A Q = Ci^Ci+x,! + Ci_ X) j ± Ci + y ; j ± Ci_y ; j), where a = s corresponds to the (+) signs and 
defines an extended s-wave operator, while a = d corresponds to the (— ) signs in the in- 
direction, and defines a d x 2_ y 2 pair operator, x, y are unit vectors along the axis, and the 
normalization is chosen such that < Z 2 h < 1- The size dependence of Z 2 h shown in Fig. 2a 
seems smooth and flat both in the large and small coupling regions suggesting that Z 2 h is 
nonzero in the bulk limit for all finite values of the coupling J/t. Fig. 2b illustrates the 
coupling dependence of Z 2 h. The results for different cluster sizes approximately follow a 
linear behavior, Z 2 h ~ J/t, in the interval < J/t < 1. This result is in agreement with 
calculations carried out in the one hole sector that suggested Z lh ~ (J/t) 1 / 2 (see ref. |p4fl ). 

To complete our study, let us consider the spectral decomposition of the pairing operator 
that can be carried out using standard techniques. It is defined as, 

P(u) = £ KV>f |A^K>| 2 % - {Ef - E%)) (0.3) 

n 

where the notation is standard. The pair spectral functions shown in Fig. 3 illustrate 
some of the conclusions of this paper. The calculations shown in the figure correspond to a 
\^20 x \^20 cluster at J/t = 0.3 and two different symmetry operators, namely d x 2_ y 2 and 

a 



extended-s wave. The difference between the two is clear. While the spectral decomposition 
of the d-wave operator shows a clear sharp peak at the bottom of the spectrum with an 
intensity given by the Z 2 h presented in Fig. 2, the decomposition for the s-wave shows no 
appreciable spectral weight at low frequencies. The fact that Fig. 3 reproduces most of the 
qualitative features found on smaller systems J7J gives credibility to small cluster calculations. 
The study of the size dependence of the d-wave pairing spectral function discussed before 
[ |16|1 shows that, with increasing system size, (i) the low energy QP peak survives, and (ii) 
the higher energy peaks eventually merge into a continuous background. 

Lastly we discuss the nature of the pairing in momentum space. The form of the d- 
wave pair operator in k-space, = J2ki cos ^x — cos ^y) c k,t c -k,|' suggests that the pairing 
occurs predominantly between two holes at momenta (71", 0) (or (0, 7r)) in agreement with the 
calculation of the single hole spectral function in one hole-doped clusters. 



Summarizing, in this paper we have presented a complete study of the behavior of two 
holes injected on a spin-1/2 antiferromagnetic background. In agreement with previous 
numerical studies carried out by the authors, we observed the tendency towards the formation 
of bound states for couplings larger that some finite critical value J/t| c . The bound state 
in that region has d x 2_ y 2 symmetry, spin zero and carries a nonzero overlap with the state 
obtained by applying a local pair creation operator over the ground state of zero hole. Such 
a behavior, signalled by a nonzero Z 2 h, favors the interpretation of the two holes bound state 
as a quasiparticle of charge 2e and spin zero, which would be the actual carriers of charge 
under an applied electric field. The analysis of the spectral decomposition of the pairing 
operator illustrates a dramatic difference between the favored d-wave bound state and an 
extended s-wave state. These quasiparticles are natural candidates to Bose condensate at 
low temperatures into a superconducting d-wave ground state in the t — J model, similar to 
that recently found by two of us (E.D. and J.R.) away from half-filling. || Such a condensate 
may become a realization of the recently proposed new theoretical ideas |T|] to explain the 
behavior of the actual high-Tc cuprate superconductors. 
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FIGURE CAPTIONS 

Figure la 

Average hole-hole distance in the ground state as a function of the coupling constant for the 
four clusters studied in this paper. Open triangles denote results for a cluster of 16 sites, 
open squares for 18 sites, full triangles for 20 sites, and full squares for 26 sites. 
Figure lb 

Binding energy, A#, of two holes in the t — J model as defined in the text. Open triangles 
denote results for a cluster of 16 sites, open squares for 18 sites, full triangles for 20 sites, 
and full squares for 26 sites. The points with the error bars joined by a dashed line are 



Green's Function Monte Carlo results taken from Ref. ||21|| . 
Figure 2a 

The wave-function renormalization, Z2h (defined in the text), as a function of the inverse of 
the number of sites N for several values of the coupling J/t between 5 and 0.1. 
Figure 2b 

The wave-function renormalization, Z2h (defined in the text), as a function of the coupling 
J/t. Open triangles denote results for a cluster of 16 sites, open squares for 18 sites, full 
triangles for 20 sites, and full squares for 26 sites. 
Figure 3 

The dynamical response of the pairing operator defined in the text obtained on a x \^20 
cluster, and at J/t = 0.3. In (a) the pairing operator used has d-wave symmetry, while in 
(b) it corresponds to extended s-wave. 
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